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The full relaxed-density potential between spherical nuclei is considered as a sum of the macro¬ 
scopic and shell-correction contributions. The macroscopic part of the potential is related to a 
nucleus-nucleus potential obtained in the framework of the extended Thomas-Fermi approach with 
the Skyrme and Coulomb forces and the relaxed-density ansatz for evaluation of proton and neu¬ 
tron densities of interacting nuclei. A simple prescription for the shell-correction part of the total 
potential is discussed. The parameters of the shell-correction and macroscopic parts of the relaxed- 
density potential are found by fitting the empirical barrier heights of the 89 nucleus-nucleus systems 
as well as macroscopic potentials evaluated for 1485 nucleus-nucleus systems at 12 distances around 
touching points. 

PACS numbers: 24.10.-i, 25.70.-z, 25.70.Jj 


I. INTRODUCTION 

Knowledge of the nucleus-nucleus potential around the 
barrier is necessary for the analysis of heavy-ion fusion 
and other nuclear reactions |.1-10|. Heavy-ion fusion at 
low energies is very important in the physics of stars 
@,UM2. As a result, the nucleus-nucleus interaction 
potential is a key ingredient for the description of various 
nuclear reactions in the nature. 

The microscopic evaluation of a potential between nu¬ 
clei is based on both the effective nucleon-nucleon force 
and the nucleon density distributions of interacting nu¬ 
clei; see Refs. [214-il IfiUot ITdl - [2f^ and papers cited therein. 
There are also simple phenomenological parametrizations 
of the nucleus-nucleus potential QlJMlIffl, which 
are often used for the description of various heavy-ion re¬ 
actions. The parameters of phenomenological potentials 
depend smoothly on the numbers of nucleons of inter¬ 
acting nuclei. When the nuclei approach each other, the 
energies of the nucleon single-particle levels of each nu¬ 
cleus are shifted and split due to the interaction of nu¬ 
cleons belonging to different nuclei. Therefore, the shell 
structures of both nuclei are changed at small distances 
between the nuclei. 

The nucleon single-particle motion and shell effects are 
studied in the framework of detailed microscopic theo¬ 
ries, which take into account the time-dependent trans¬ 
formation of the two-center field into a single-center one 
HH-HU and the dissipation of the collective motion. 
Such microscopic approaches are related to the complex 
numerical calculations. The goal of this study is to find 
a simple and easy-to-use approximation for the nucleus- 
nucleus potential, which takes into account the bulk and 
shell effects. 

The full energy of a nucleus consists of the sum of 
macroscopic and microscopic contributions according to 
the shell-correction method proposed by Strutinsky [37|. 
The shell effect contribution has been successfully eval¬ 
uated using the Strutinsky shell-correction prescription 
[37l HH, which is simpler than the full microscopic treat¬ 


ment. The shell-correction energy is an important in¬ 
gredient of various modern nuclear physics models be¬ 
cause the shell corrections strongly improve the modeling 
accuracy. Many key parameters such as atomic masses 
[39ll40l fission process [UHl] , cluster decay of heavy nu¬ 
clei [u], |42| , and fusion cross section hindrance at deep 
sub-barrier energies [43) are well described in the frame¬ 
works of the macroscopic-microscopic models. 

The contribution of the shell structure to the nucleus- 
nucleus potential has fully been ignored in phenomeno¬ 
logical approaches. Therefore, it is desirable to find an 
improved phenomenological potential which takes into 
account the shell structure variation induced by the inter¬ 
action between nucleons in the interacting nuclei. Such 
a potential should be more precise because it takes into 
account the unique shell-correction contribution of inter¬ 
acting nuclei. The shell-correction energy varies from 
one nucleus to another substantially [371 ). In contrast, 
the macroscopic part of the full potential is described 
by a simple phenomenological expression which depends 
smoothly on the numbers of protons and neutrons of in¬ 
teracting nuclei. So, the macroscopic-microscopic poten¬ 
tial should take into account both gross and individual 
properties of the specific nucleus-nucleus system. 

A simple approach to the full nucleus-nucleus poten¬ 
tial consisting of macroscopic and microscopic parts for 
a wide range of nucleus-nucleus systems is developed in 
the present paper. The parameters of the full potential 
are found using both data for empirical nucleus-nucleus 
barriers and the values of macroscopic nucleus-nucleus 
potentials around the barrier for various systems. The 
compilation of the data for empirical nucleus-nucleus bar¬ 
rier heights is given in Refs. |36l . l44l - l47| . The macroscopic 
nucleus-nucleus potential is evaluated in the framework 
of the non-local extended Thomas-Fermi (ETF) energy- 
density functional theory and the relaxed-density ap¬ 
proach for the description of the nucleon density distri¬ 
bution of two interacting nuclei around the touching dis¬ 
tance. The energy-density functional takes into account 
the contributions of the nucleon kinetic energy and the ef- 
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fective nucleon-nucleon Skyrme and Coulomb forces |48| . 
The ETF kinetic energy-density functional consists of the 
Thomas-Fermi term and all correction terms of the or¬ 
der ft 4 (48|- Such an approach shows high accuracy for 
the description of atomic masses 39, 48], nucleus-nucleus 
potentials for various systems j 6 i; , | 20 |, | 22 ( and the cluster 
decay of heavy nuclei 421 ■ 

A detailed description of the macroscopic and shell- 
correction contributions to the full nucleus-nucleus po¬ 
tential is given in Sec. 2. The shape and parameters 
of the relaxed-density nucleus-nucleus potential with the 
shell-correction contribution are defined in Sec. 3. Sec. 4 
is devoted to a discussion of the results and conclusions. 


II. RELAXED-DENSITY NUCLEUS-NUCLEUS 
POTENTIAL 


In the shell-correction method ( 37 ] , the full interaction 
potential energy of the system of interactin g nu clei is 
written in the form (see also Refs. ful Id 3 l l 49 tl 50 |l 

Hot (R) = V macro (R) + V sh (R) (1) 

= [E\ 2 (R) — E\ — E 2 ] 

+ [ 6 E 12 {R) - 5E! - 5E 2 }. 


Here E\, E 2 , 5E\ and 8 E 2 are the macroscopic and shell- 
correction energies of the non-interacting spherical nuclei 
1 and 2 , respectively. E\ 2 {R) and 5Ei 2 {R) are the macro¬ 
scopic and shell-correction energies of the interacting nu¬ 
clei at distance R between the mass centers of separated 
nuclei. The shell-correction energies 5 E\ 2 {R), 5E\ and 
8 E 2 include proton and neutron shell correction energies 
related to both the non-uniformity of single-particle spec¬ 
tra around the Fermi energies and the pairing corrections 


A. Macroscopic interaction potential between 
nuclei 

The macroscopic part of the interaction potential en¬ 
ergy Hnacro (R) between spherical nuclei is given by 

H macro (i?) = E 12 (R) -E 1 -E 2 . ( 2 ) 

The macroscopic energies of nuclei in Eqs. (1) and (2) are 
evaluated with the help of a semiclassical ETF energy- 
density functional £ [p p (r), p„(r)]. These energies are 

Eu(R) = J £[p p (r,R),p n (r,R)]dr, ( 3 ) 

E k = j £[pkp(r),pkn(r)\ dr, ( 4 ) 

where p\ p (r), p 2 P (r), pi n {r) and p2n{r) are the proton 
and neutron densities of the non-interacting nuclei, while 
p p (r,R) and p n (r,R) are the proton and neutron densi¬ 
ties of the interacting nuclei, k = 1 , 2 . 


1. Energy-density functional 


The semiclassical energy-density functional includes 
the Skyrme and Coulomb interactions as well as the in¬ 
trinsic proton and neutron kinetic energies obtained in 
the ETF approach [48|. According to J4|| the following 
expression for the energy density functional has been de¬ 
duced: 

£[Pp(r),p n { r)] = ^-[r p (r) +r n (r)] + Vsk(r) + V c (r).( 5 ) 

The kinetic parts for protons (* = p) and neutrons (i = n) 
are given by 


n( r) = T iT F{r) + n 2 (r) + 774(1-), (6) 


where TiTp{r) is the Thomas-Fermi contribution to the 
kinetic energy-density functional, and 772 ( 1 ") and 774 ( 1 ") 
are semiclassical ti 2 and H A correction terms to the kinetic 
energy-density functional for the non-local case, respec¬ 
tively. 

The nuclear interaction part Vsk(r) results from the 
Skyrme force and reads 

Vsk(r) = y ]{l + ^x 0 )p 2 - {x 0 + ^){p 2 p +p 2 n )\ (7) 

+ J^ t 3P °‘[( 1 + \x 3 )p 2 - (x 3 + ^)(pp + pl)\ 

+ jIM 1 + \ x i) + ^(1 + ^x 2 )]rp 

+ ^[h{x 2 + 7 ,) ~ tl(X! + ^)}(t p Pp + T nPn ) 

+ ^[ 3 H(i + ^l) - t 2 ( 1 + \x 2 )](ypf 

~ + ~j) + ^( X 2 + ~ )]((Vp„) 2 + (V/3p) 2 ) 

Wn 

-y [J ' Vp + J„ ■ Vp n + J p ■ Vp p ] , 


where to, t\, t 2 , Xq, X\, X 2 , a and Wq are the Skyrme-force 
parameters, J; are the spin-orbit densities, p = p p + p n , 
t = Tp + T n and J = J p + 3 n . The expressions for 77 2 ( r) , 
774 ( 1 -), and Jb for a non-local case are given in Ref. 481. 

The Coulomb-energy density is given by 


V C (r) = — Pp (r) 


/ 


P P (r') 


dr' 




1/3 

M r )] 4/3 , 


( 8 ) 


where the first term is the direct contribution, the second 
term is the local approximation to the exchange contri¬ 
bution, and e is the proton charge. 

This approximation for the energy of nucleus based on 
the semiclassical energy-density functional is simple and 
accurate. Similar approximations ware successfully used 
for the evaluation of the atomic masses [3!|, the fission 
barrier characteristics [48| , the nucleus-nucleus pot entials 
a da da sa mm and the cluster emission [42j ■ 
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2. Parametrization of density distribution 

The nucleon density is a necessary ingredient of the 
semiclassical energy-density approach; see Eqs. (3)-(8). 
It is difficult to find the proton and neutron density dis¬ 
tributions by solving the integro-differential variational 
Lagrange equations in the framework of the non-local ft 2 
ETF approach for the case of a spherical nucleus 33- 
The neutron and proton densities of nuclei in the case of 
the non-local h A ETF approach are found by using trial 
functions only [48;]. The proton and neutron density dis¬ 
tributions for system of interacting nuclei have not been 
evaluated yet in the framework of the ETF. 

As a rule the densities of interacting nuclei are 
parametrized according to specific physical conditions of 
the reaction. The frozen approximation for the density 
distributions of interacting nuclei is often used for the 
evaluation of the nucleus-nucleus potentials in the frame¬ 
works of the energy-density IE HE HE HE . I2E I2E] and 
double-folding Hi M HE Hlj approaches. Note that 
it is reasonable to use the frozen approximation for fast 
nucleus-nucleus collisions, when the proton and neutron 
densities cannot significantly relax during the collision 
process @. The proton or neutron densities at the fixed 
point of space in this approximation are the sum of the 
corresponding nucleon densities of each nucleus in this 
point. As the result, the nucleon density can exceed the 
equilibrium density of the nuclear matter p in some space 
region at small distances between nuclei. 

The nucleus-nucleus fusion reactions around the bar¬ 
rier and, especially, deeply below it are related to the slow 
process of collective motion of nuclei. During this pro¬ 
cess the proton and neutron densities can relax and the 
densities of nuclei cannot be simply presented as sums of 
nucleon densities of the interacting nuclei. The relaxed- 
density distributions should fulfill obvious conditions: 

- The values of density at any point in space do not 
exceed the equilibrium density of nuclear matter 
p because the compressibility modulus of nuclear 
matter strongly prevents excess of p. 

- The values of relaxed density at any point in space 
should be smaller than the one in the frozen ap¬ 
proximation but larger than the density values of 
any of the interacting nuclei at this point. 

Taking into account these conditions, the relaxed-density 
ansatz [42j for the proton (neutron) density of the inter¬ 
acting nuclei at point r can be presented as 

Pp(n)(p,R) = Pip(„)(r) + P2p(n) (E R) (9) 

2 Plp(n)(r) P2p(n) (E R) 

Plp(n ) T P2p(n) 

Here Pip(n)(^) = Pip(n)fip(n)(p') and p 2 p(n)(E-^) 

P 2 p(n)f 2 p(n) (e R) are the density distributions of the 
proton (neutron) subsystem related to the first and 
second nuclei, respectively. The functions /i p ( n )(r) 


and f 2p { n ) ( r , R) describe the space distributions of cor¬ 
responding densities. These functions are normal¬ 
ized to unity at the center of corresponding nuclei 

/lp(n) (**) |r=0 / 2 p(n) (E R) |r=0 L 

Let us consider the cases of strongly and weakly over¬ 
lapping nuclei in the framework of parametrization (9). 
For the sake of simplicity, similar shapes of the proton 
and neutron density distributions in each nucleus are pro¬ 
posed for consideration of these cases. Therefore, the 
density distributions are 

PipO) = ^P/iO), Pin(r) = ^p/i(r), 

Zo No 

p 2 p (r, R) = — p/ 2 (r, R ), p 2 n (r, R) = — 1 p/ 2 (r, R), 

>12 Vi 2 

where Zi and N t are the numbers of proton and neutron 
in the nucleus i, Ai = Zi + Ni , i = 1 , 2 . 

The nuclei overlap strongly at small distances between 
them. The total densities of each nucleus at the point 
close to the center of each nucleus are 

Pi(M) =p(i-si), 
p 2 (|r|, i?) = p 2 (R- |r|) = p(l - 52 ), 

where 51 = 1 - /i(|r|) < 1 and g 2 = 1 - f 2 (\r\,R) < 1. 
Substituting these densities into (9) I find the value of 
relaxed density at small distance between nuclei at the 
point between the centers: 

P (E R) ~ P[ 1 - (fli + 92 )/2] < p. 

In contrast to this, the total density of these nuclei at the 
same point in the frozen density approach is 

p froZ en(r,i?) = Pi(|r|) + p 2 (|r|,-ff) =p[ 2 - gi - g 2 \ >p. 

Note that 

p{ r, R) « p[ 1 - (31 + <? 2 )/2] < Pfrozen(E R) ■ 

In the opposite case of well separated nuclei (large dis¬ 
tance between interacting nuclei) the nucleon densities of 
each nucleus in the point between nuclei and far from the 
center of each nucleus are very small 

Pi(M) = P/i(l r l) = Pfi, 

p 2 (|r|,i?) = p/ 2 (|r|, R) =pf 2 , 

where /^ 1. Substituting these densities into (9) I 

obtain the value of relaxed density at this point 

p(r, R) « p(fi + f° ~ 

Comparing this density with the other ones, it is easy to 
conclude that 

- The value of relaxed density is larger than the den¬ 
sity of each nucleus at this point, i.e. 

p(r,R) > pi(|r|) =pfi, p(r,R) > p 2 (|r|,i?) =pf 2 ■ 
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The value of relaxed density is smaller than the sum 
of densities of these nuclei at this point, i.e. 

p( r, R) ~ p(fi + /° - /i 0 / 2 °) < Pi(r) + P 2 (r, i?) 

= P(/l + / 2 °) = Pfrozen(r, A). 

Thus, the parametrization (9) satisfies the proposed 
conditions. This parametrization leads to a more realis¬ 
tic spatial distribution of density formed by two nuclei 
near the touching point at collision energies around the 
barrier height than the one obtained in the frozen density 
approach. Approximation (9) also drastically simplifies 
the numerical calculations of the macroscopic potential 
Vm (R) in the framework of the ETF approach with the 
Skyrme and Coulomb forces. 

The proton (neutron) density distributions of sepa¬ 
rated spherical nuclei are evaluated in the framework of 
the Hartree-Fock method using the SLy4 parameter set of 
the Skyrme force [541. The BCS pairing with the Lipkin- 
Nogami method [531 ] for approximate particle number 
projection is used in the calculation. The microscopic 
proton and neutron densities are given on the mesh with 
step h. At small distances between nuclei, when the den¬ 
sities of nuclei are overlapped, let us propose that the 
proton and neutron density values of each nucleus are 
the same as before at any mesh points. However, the 
size of the mesh is slightly changed to a p ( n )h. The scal¬ 
ing factors a p („) can be fixed at any value of R by the 
conservation conditions of the total numbers of protons 
Z 1 +Z 2 and neutrons N 1 +N 2 in the system of interacting 
nuclei, which are 

Jdr p p („)(r, R) = Z\ + ^ 2 (= Ni + -/V 2 ). (10) 

The values of a p ( n ) are very close to 1 at R around the 
touching point of the nuclei. 

By using ansatz (9) for the relaxed-density 
parametrization, the microscopic Hartree-Fock den¬ 
sities and conditions (10), and the description of 
densities of two nuclei at large R and around the 
touching point is simplified. As a result, the macroscopic 
interaction potential energy of the nuclei 14iacro(l?) can 
be easily obtained by substituting these densities into 
Eqs. (2)-(8). 

B. Shell-correction contribution to the full 
interaction potential between nuclei 

The shell-correction contribution to the full interaction 
potential energy is 

V sh (R) = 6E 12 (R) - 5E 1 - 6E 2 . (11) 

It is obvious that mutual influence of nuclei on their 
single-particle spectra is negligible at distances much 
greater than the sum of radii of nuclei Rt = R\ + i? 2 , 
therefore 

( 12 ) 


and 

V s h(R)\R»Rt = 0 . 

The SE 12 ( R ) value at R « 0 equals the shell-correction 
energy of the nucleus formed in the fusion of nuclei 1 
and 2 at the equilibrium deformation. The full potential 
should be equal to the fusion reaction Q-value at R = 0, 
i.e. Vtot(R = 0) = Q. 

The nuclei strongly interact at small distances between 
them. This interaction leads to the shift and the split¬ 
ting of the single-particle levels in both nuclei. Due to 
this the proton and neutron single-particle spectra in the 
vicinity of the respective Fermi levels become more ho¬ 
mogenous around touching distance Rt- Such behavior 
of the single-particle levels is clearly demonstrated in the 
frameworks of the two-center shell model [ 3 , [§, HH, US] 
and microscopic cluster decay model |4l[ . 

According to the shell-correction method [13, SI], the 
absolute value of the shell correction energy is reduced 
in the case of the more homogeneous single-particle spec¬ 
trum in the vicinity of the Fermi level. The sharp re¬ 
duction of the shell-correction contribution to the full 
potential energy around the touching point of nuclei is 
obtained in the framework of the microscopic cluster de¬ 
cay model [4l| . So, the perturbation of the single-particle 
nucleon levels, which is proportional to the strength of 
nucleon-nucleus interaction, reduces the absolute value 
of the shell-correction energy. 

The perturbation of the single-particle level energies 
and the interaction between nucleons belonging to dif¬ 
ferent nuclei are increased when the distance between 
the surfaces of the interacting nuclei decreases. The per¬ 
turbation strength is related to the density distribution 
in the nucleus induced by the disturbance as well as to 
the radius of the nucleon-nucleon force. The density dis¬ 
tribution is often parametrized by the Fermi distribu¬ 
tion. Since the radius of the nucleon-nucleon force is very 
short, it is reasonable to approximate the shell-correction 
contribution to the full nucleus-nucleus potential around 
the touching point of two nuclei as [Hj, [43| 

5E 12 (R) « [SEx + SE 2 \f sh {R). (13) 

Here 

f sh (R) = 1/ {1 + exp [(R sh - R)/d sh ]} , (14) 

and R s h and d s h are the radius and the diffuseness re¬ 
lated to the attenuation of the shell-correction magnitude 
with reduction of distance R. This approximation for the 
shell-correction energy contribution into the full potential 
is rough, but it can greatly simplify the calculations of 
the shell-correction energies for various nucleus-nucleus 
systems at different distances R. Note that numerical 
calculations of the shell-correction energy in the frame¬ 
work of the two-center model (4lJ show that the values 
of the shell-correction energy of the nucleus-nucleus sys¬ 
tem abruptly change around the touching point of two 


<5Ei 2 (A)|_R>_R t = 8 E\ + 8E2 


5 


nuclei. Therefore, the radial dependence of the shell- 
correction energy described by Eqs. (13) and (14) is 
reasonable. Moreover, the exponential reduction of the 
shell-correction energy values related to washing out the 
non-homogeneity of the single-particle spectra is often 
considered in nuclear physics 0, S Mm- 

The shell-correction energy of interacting nuclei 
6 Ei 2 (R) smoothly approaches the limit of non¬ 
interacting nuclei (12) at large distances R between nu¬ 
clei, i.e., V s h(R) —> 0 at R R s h ~ Rt- As a result, the 
shell-correction contribution to the full nucleus-nucleus 
potential is 


V sh (R) - [SE 1 + SE 2 } • [f sh {R) - 1], (15) 

This approximation for the shell-correction contribution 
to the potential energy of two nuclei is successfully used 
for the description of both the cluster emission from 
heavy nuclei [421 ] and the fusion cross section hindrance 
at deep sub-barrier energies [43| . 

It is possible to evaluate the shell-correction energy val¬ 
ues 5Ei and 6 E 2 according to the Strutinsky prescription 
[371138 ] for a specific nucleon mean-field. However, the 
easiest way to estimate the value of the shell-correction 
energy SE in a spherical nucleus is to find the differ¬ 
ence between the experimental B exp and macroscopic B m 
binding energies of the nucleus: 

6 E = B exp - B m . (16) 


Moreover, this approach for the evaluation of the shell- 
correction energy is similar to the one applied for the 
evaluation of the energy level density in nuclei; see Refs. 

[H®. 

The values of B exp in Eq. (16) can be found in the re¬ 
cent evaluation of the atomic masses [Hj], while the value 
of B m in the nucleus with Z protons and N neutrons is 


B m = -15.86864A+ 21.18164A 2/3 - 6.49923A 1/3 + 


'N — Z 

2 

26.37269A - 23.80118A 2/3 - 8.62322A 1/3 

A 


J 


+ 


Z 2 
AV3 L 


0.78068 — 0.63678A -1/3 


+ P p + P n . (17) 


Here B m is the binding energy in MeV, A = Z + N, P p ( n ) 
are the proton (neutron) pairing terms, which equal to 
P p(n) = 5.62922(4.99342)A -1 / 3 in the case of odd Z (N) 
and Pp( n ) = 0 in the case of even Z ( N ). Equation (17) is 
a simple extension of the liquid drop Weizsacker formula 
|62l . 1631 ] for the binding energy of nuclei. The coefficients 
in Eq. (17) are evaluated by fitting the recent values of 
the atomic masses [61]. The experimental binding ener¬ 
gies of 3353 nuclei are described by Eq. (17) with a root 
mean error of 2.49 MeV. This error is very small com¬ 
pared to the experimental values of the nuclear binding 
energies in medium and heavy nuclei. 

The liquid drop mass formula may be found by us¬ 
ing the Skyrme force too [64 ]. Such an approach may 


look more consistent because the macroscopic nucleus- 
nucleus interaction energy is evaluated with the help of 
the Skyrme energy density functional; see Eqs. (2)-(8). 
However, the liquid drop mass formula obtained in Ref. 
[64| using the Hartree-Fock calculation with the SLy4 
parametrization of the Skyrme force leads to a larger 
value of the root mean error for experimental masses 
than the one in Eq. (17). The atomic masses of very 
heavy nuclei calculated in the framework of the Hartree- 
Fock approach with the SLy4 parameter set are system¬ 
atically underbound Therefore, Eqs. (16) and (17) 
are applied for the evaluation of the shell-correction en¬ 
ergy here because it is better to use a simple and high 
accuracy approach. 

Note that slightly different values of the shell- 
correction energy for the same nucleus can be calculated 
applying various formulas for the liquid drop mass for¬ 
mula [601 ] or using diverse parametrizations of the nu¬ 
cleon mean field for the shell-correction calculations ac¬ 
cording to the Strutinsky prescription. Similarly, the 
Hartree-Fock calculations with various parameter sets of 
the Skyrme force lead to different values of the shell- 
correction energy for the same nucleus. For example, 
the values of the shell-correction energy of 208 Pb eval¬ 
uated by using Eqs. (16)-(17), the Megnoni-Nakamura 
approach [60t l. and the Nilsson and Saxon-Woods mean- 
field models [38l. l66j] are —10.620, —9.972, —12 and —12 
MeV, respectively. These different values of the shell- 
correction energy are slightly affected by the value of the 
nucleus-nucleus potential due to the radial dependence of 
the shell-correction contribution, see Eqs. (14) and (15). 
For example, the variations of the nucleus-nucleus po¬ 
tentials for the system 16 O+ 208 Pb at the touching point 
induced by various values of the shell-correction energies 
of 208 Pb (-10.620 and -9.972 MeV or -10.620 and 12 
MeV) are 0.16 or 0.34 MeV, respectively. Moreover, the 
corresponding differences of the nucleus-nucleus poten¬ 
tials near the barrier distances are even smaller than 0.16 
or 0.34 MeV due to the exponential dependence f s h{R) 
on R. Note that the value of the barrier height of the 
nucleus-nucleus potential for the system 16 O+ 208 Pb is 
close to 77 MeV [43j. Therefore, the mo del-dependent 
variation of the shell-correction energy is not crucial for 
the accuracy of the nucleus-nucleus potential. 

Any influence of the relative motion of the nuclei on the 
nucleus-nucleus interaction and the single-particle levels 
has been ignored up to now. However, the nucleons move 
in the approaching nuclei during reaction and the interac¬ 
tion of nucleons belonging to different nuclei disturbs the 
nucleons in both nuclei. This leads to time-dependent 
disturbance of the shell structure of each nucleus. It 
is obvious that the strength of this disturbance should 
depend on the ratio of the nucleon velocity (kinetic en¬ 
ergy) and the relative velocity (kinetic energy) of the nu¬ 
clei. The relative velocity of two nuclei depends on ra¬ 
dial and rotational degrees of freedom. When the relative 
velocity of approaching nuclei is small, the static (adia¬ 
batic) consideration is close to realistic because the shell 
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structure of the nuclei can be relaxed due to high veloc¬ 
ities of the nucleons in the nuclei. In contrast, the shell 
structure and nucleon density distributions cannot be dis¬ 
turbed at fast collisions. Therefore, nuclei can touch each 
other without modification of the shell structure of col¬ 
liding nuclei at high collision energies. This conclusion 
is supported by time-dependent microscopic calculations 
pil [ 25 I ] , which show that the nuclei overcome the barrier 
essentially in their ground-state density at high collision 
energies. Thus, the shell-correction contribution to the 
full potential depends on the collision energy E. The en¬ 
ergy dependence of this contribution can be presented in 
the form [43[ 


V sh (R, E ) 


V sU R ) -exp[-a{E-B)}, at E > B,. . 

V° h (R), at E < By*’ 


where B is the barrier height of the macroscopic poten¬ 
tial Vtnacro(-R) and a is the reduction parameter of the 
shell-correction energy contribution to the full potential. 
So, the shell-correction contribution to the full potential 
decreases with increasing the collision energy. Eq. (18) 
does not depend on the velocity of the nucleons because 
nucleon velocity at the Fermi level has a small variation 
from one nucleus to another and the value of the shell- 
correction energy depends on the inhomogeneity of the 
single-particle spectra around the Fermi energy [37j |. 

The full potential at E > B moves smoothly to¬ 
wards the macroscopic potential with increasing E. The 
nucleus-nucleus collision process is near adiabatical at 
small collision energies E < B, because the velocity of 
relative motion is less than the velocities of nucleons in 
nuclei. Therefore, the energy dependence of the shell- 
correction contribution is ignored at low collision ener¬ 
gies. 

The value of the parameter a = 0.35 MeV -1 in Eq. 
(18) is determined by fitting the fusion cross section for 
the reaction 16 0 + 208 Pb [43|] because it has been mea¬ 
sured in a very wide energy range around the barrier 
[67}. The thicknesses of the barrier at deep sub-barrier 
collision energies become larger due to the positive shell- 
correction contribution to the full potential for the re¬ 
action 16 0 + 208 Pb. This directly leads to a reduction 
of the barrier penetrability and, as a result, to the hin¬ 
drance of the deep sub-barrier fusion which is observed 
for this reaction 0, iTfT . 


III. PHENOMENOLOGICAL 
RELAXED-DENSITY NUCLEUS-NUCLEUS 
POTENTIAL 


The proton and neutron densities of the nuclei are ob¬ 
tained in the framework of the Hartree-Fock model and 
the relaxed-density ansatz (9). The macroscopic relaxed- 
density nucleus-nucleus potential is calculated substitut¬ 
ing these densities into Eqs. (2)-(8). The parameter set 
SLy4 [54J of the Skyrme force is used for the evaluation of 


both the Hartree-Fock nucleon densities and the macro¬ 
scopic part of the nucleus-nucleus potential. This param¬ 
eter set is adopted in the calculations because the nuclei 
along and far away from the beta-stability line are often 
used in heavy-ion reactions and this set of the Skyrme 
force is adjusted to reproduce various properties of nu¬ 
clear matter, neutron stars, and finite nuclei along and 
far away from the beta-stability line. 


The values of fusion barrier heights obtained in 
the macroscopic relaxed-density approach for light and 
medium nucleus-nucleus systems are close to the empiri¬ 
cal ones from Refs. J3g,|45|, |47|. Unfortunately, the eval¬ 
uated values of the fusion barrier heights for heavy and 
very heavy nucleus-nucleus systems are higher than the 
empirical ones. Note that a similar tendency is found in 
Ref- 01 w h en the barrier heights of nucleus-nucleus po¬ 
tentials are obtained in the framework of the semiclassical 
energy-density approach with the frozen density distribu¬ 
tion of two nuclei [35| and the SkM* parameter set of the 
Skyrme force @ (see also Fig. 1). Such a tendency is 
probably related to a systematical underbounding |65( of 
very heavy nuclear systems and/or nuclei in the frame¬ 
work of the Skyrme-Hartree-Fock approach. Moreover, 
the overestimation of empirical barrier heights for heavy 
and very heavy nucleus-nucleus systems is observed for 
various phenomenological potentials too; see Fig. 1. 


Our goal is to find the phenomenological relaxed- 
density potential which describes well the empirical bar¬ 
riers for light, medium, heavy, and very heavy nucleus- 
nucleus systems. The values of the macroscopic part of 
this potential around the touching distances should be 
close to the numerically evaluated relaxed-density macro¬ 
scopic potential. To solve this task the next approach is 
applied. 


At the beginning, 54 spherical nuclei between 16 O 
and 216 Po are selected. Then 1485 macroscopic relaxed- 
density potentials between all possible pairs of these 54 
nuclei at 12 distances around the touching point are eval¬ 
uated. These distances equal Rij = 1.18(Aj{ 3 + A \[ 3 ) + 
0.3j—0.6 frn, where An and A. 2 i are the numbers of nucle¬ 
ons in the interacting nuclei, * = 1,2, 3,..., 1485 and j = 
1,2,3,..., 12. The dataset for 17820 (17820 = 1485 x'l2) 
macroscopic potential values Vi maclo (Rij) and the data 
for 89 empirical barrier heights B^ mp between spheri¬ 
cal nucleus-nucleus systems from Refs. 36f, HH-Hg are 
used for the evaluation of the macroscopical and shell- 
correction parts of the full phenomenological relaxed- 
density nucleus-nucleus potential. Note that, as a rule, 
the values of the empirical barrier heights are extracted 
from the analysis of subbarrier heavy-ion fusion reac¬ 
tions. Therefore, the shell-correction contributions into 
the full nucleus-nucleus potentials are evaluated in the 
case E < B. The parameters of the phenomenological 
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potential are found by the minimization of the function 


F 


15 


1 

17820 


1485 12 


EE 


(yphen , R s 
[ v ^macro'- > 

y ^macro(^j) 



89 

Y.{BT - sf")' 


fc=l 


1/2 


21 V 2 


(19) 


Here R£ hen are the values of the barrier heights evalu¬ 
ated by using the full phenomenological potential. The 
values of the macroscopical part of the phenomenological 
potential should be similar to the ones of the numerically 
evaluated relaxed-density potential due to the minimiza¬ 
tion of the first term in Eq. (19). The minimization of 
the second term of Eq. (19) leads to correct values of 
barriers for various nucleus-nucleus systems. The factor 
15 in the first term of Eq. (19) enhances the influence 
of fitting of the macroscopic potentials, and the values of 
the first and second terms of Eq. (19) are close to each 
other in this case. 

The full phenomenological relaxed-density nucleus- 
nucleus potential is the sum of the Coulomb and nuclear 
parts 


V tot (R) = Vc(R) + Enuci (R), 


( 20 ) 


where 

V C (R) 


( Z 1 Z 2 e 2 


1 R ’ 



3 in 2 

l Rc 

L 2 2 


at R > Rc, 
at 0 < R < Rc, 


( 21 ) 


Kuci (R) 


^macro(^) + bsh(-R)j &t R > Rt, 
V™ cl (R), at0<i?<i? t , W 


and 

Ri = 1.2536H- /3 - 0.80012H“ 1/3 - 0.0021444/H, (23) 


are the radii of corresponding nuclei in frn and i = 1,2. 
The expression for the Coulomb potential is standard 
for the theory of heavy-ion reactions M- Let us put 
Rc = Rt = R\ + 7?2 for the sake of the reduction of the 
parameter number. 

The macroscopic nuclear part l^(ac ro (7?) of the poten¬ 
tial is taken in the Woods-Saxon form 


KJacro (R) 


Vo 

1 + exp[(i? - R t )/d] ’ 


(24) 


where 


Vo = Vl C + v 2 C 1/2 (25) 

is the potential strength in MeV, v\ = —27.190 MeV 
fm" 1 , Vl = -0.93009 MeV fm" 1 / 2 , 


(26) 


is the diffuseness in fm, d± = 0.78122 fm, d 2 = —0.20535 
fm 2 and C = R 1 R 2 /Rt is in fm. 

The shell-correction contribution V s h(R) to the total 
potential at R > R t is evaluated using Eqs. (14)-(18). 
The shell-correction contribution to the potential is small 
near the barrier radius and important close to the touch¬ 
ing point. Therefore, it is impossible to fix the param¬ 
eters of the shell-correction potential unambiguously by 
using our procedure. However, the values of the deep 
sub-barrier fusion cross sections are sensitive to the shell- 
correction contribution of the heavy-ion potential [43j . 
Therefore, the radius R s h and diffuseness d s h parameters 
of the shell-correction contribution [see Eq. (14)] can 
be determined using the data for the fusion cross section 
hindrance at deep sub-barrier energies for three reactions 
discussed in Ref. [HJ. The values of these parameters are 
located in narrow ranges for various fusion reactions |43| . 
So, I can use the averaging of the corresponding values 
obtained in Ref. [43;] for fixing values of R s h and d s h- The 
average value of the diffuseness parameter is d s h = 0.233 
fm. Parameter R s h can be related to the touching dis¬ 
tance of nuclei R tl i.e., R s h = R t — 0.26 fm. 

The nuclear part of the phenomenological potential at 
small distances is also taken in the Woods-Saxon form 


Kucl(^) 


Qeff 

1 + exp[(i? — i?in)/d in ] ’ 


where 


Qeff — Q 


3ZiZ 2 e 2 

2R t 


Rin ~ Rt dinlog (w) O’ 

= K 2 ucl(^t) ( Qeff _ \ 
QeffL^uc^^t) V^nucl(-Rt) / 


(27) 


(28) 

(29) 

(30) 


Here Q is the Q-value of the heavy-ion fusion reaction 
evaluated by using the recent values of the atomic masses 
[(Hj . If the binding energy of a nucleus is not given in 
Ref. [6l|, then it can be obtained with the help of Eq. 
(17). The radius R, n and the diffuseness di n of the inner 
nuclear potential are obtained by using the continuity 
conditions of both the potential 


Kucl(i?t) = VZlro(R t) + V sh (Rt) = V™ cl (R t ) (31) 


and the potential derivative 

Kuci(^t) = V^'(Rt) (32) 


at the matching point R = Rt. 

The nucleus-nucleus interaction at large distances R > 
Rt contains the macroscopic nuclear, Coulomb, and shell- 
correction contributions, see Eqs. (20)-(22). In the gen¬ 
eral case, the heavy-ion potential energy at distances 
much smaller than R t depends on the shape of fusing 
nuclei. For the sake of simplicity, the potentials in Eqs. 
(24) and (27) are parametrized by the Woods-Saxon one, 


d — d\ + d 2 /C 




















FIG. 1: ’’(Color online) Absolute differences between 89 em¬ 
pirical barrier heights [3(| I45l447| and barrier heights evalu¬ 
ated using the relaxed density with the shell-correction contri¬ 
bution (relax), frozen-density (frozen) [35), Krappe-Nix-Sirk 
(KNS) [33j , Bass (Bass80) [l}, Bass (Bass74) [27}, proxim¬ 
ity (Prox77) [30| . proximity (Prox2000l [dill . Winter [331], and 
Siwek-Wilczynska-Wilczynski (SWW) [Toj] potentials. 


because it is very often used in the theory of nuclear reac¬ 
tions. Therefore, this potential can be easily integrated 
into existing nuclear reaction codes. 

According to Eq. (1) the potential energy of a nucleus 
formed at complete fusion of two nuclei equals the re¬ 
action Q-value, i.e., Vtot(O) = Q. The relaxed density 
potential is approximately equal to Q at R = 0 due to 
Eqs. (20)-(22), and (27)-(30). 


IV. RESULTS AND DISCUSSION 

The values of 89 empirical barrier heights between 
spherical nuclei [36j, 1451 - 147) have been thoroughly de¬ 
scribed by a phenomenological relaxed-density poten¬ 
tial with shell-correction contributions. The root mean 
square error of the barrier height description related to 


FIG. 2: (Color online) Relative differences between 89 em¬ 
pirical barrier heights and barrier heights evaluated with the 
help of various potentials. The notations are the same as in 
Fig. 1. 


the second term of Eq. (19) is 0.879 MeV. The value of 
the same quantity obtained without the shell-correction 
contribution to the potential is 0.976 MeV. The compar¬ 
ison of the root mean square errors shows that the shell- 
correction contribution to the nucleus-nucleus potential 
is an important ingredient which improves the descrip¬ 
tion of empirical barrier heights. 

The root mean square errors related to the descrip¬ 
tion of 89 empirical barrier hei ghts obtained in the 
frameworks of the frozen-density 1351]. Krappe-Nix-Sirk 
311, Bass 1974 [27l . Bass 1980 [T|, proximity 1977 j3Cj, 


proximity 2000 [34j], Winter [33j, and Siwek-Wilczynska- 
Wilczynski [TO) potentials are 3.373, 3.460, 6.718, 1.815, 
1.784, 3.710, 1.835, and 1.521 MeV, respectively. These 
values of the root mean square errors are larger than the 
one obtained using the phenomenological relaxed-density 
potential with shell-correction contributions. 

The absolute differences between 89 empirical barrier 
heights and the barrier heights evaluated using various 
phenomenological potentials are presented in Fig. 1. 
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FIG. 3: (Color online) Nucleus-nucleus potentials for the sys¬ 
tem 58 Ni+ 58 Ni evaluated with the help of various approaches. 
The empirical value of the barrier height is shown by a hor¬ 
izontal arrow and the touching distance of these nuclei by a 
vertical arrow. 


As a rule, various empirical potentials well describe the 
values of the barrier heights for light nucleus-nucleus 
systems. The values of the barrier heights evaluated 
with the help of phenomenological potentials from Refs. 
Mm [33l - l35| for heavy and especially very heavy 
nucleus-nucleus systems are larger than the empirical 
ones. It is also useful to compare the relative differ¬ 
ences between the empirical barrier heights and the bar¬ 
rier heights obtained using various phenomenological po¬ 
tentials. The corresponding differences are presented in 
Fig. 2. The relative differences between the empirical 
and theoretical barrier heights clearly show the accuracy 
of the barrier evaluation for light systems in the frame¬ 
work of various approaches. 

The analysis of the results presented in Figs. 1 and 
2 shows that phenomenological relaxed-density potential 
with shell-correction contributions successfully describes 
the empirical barrier heights in the full range of nucleus- 
nucleus systems in contrast to other potentials. 



FIG. 4: (Color online) Full relaxed-density potential with the 
shell-correction contribution for the system o8 Ni+ j8 Ni evalu¬ 
ated for various contributions of the shell-correction energies 
to the potential. 


The nucleus-nucleus potentials for the system 
58 Ni+ 58 Ni evaluated in various approaches are presented 
in Fig. 3. The value of the shell-correction energy 
for 58 Ni is —4.45 MeV. The barrier height of the 
phenomenological relaxed-density potential with the 
shell-correction contribution is close to the empirical 
value 97.90 MeV [47]; see Fig. 3. The relaxed-density 
potential deviates from other potentials around the 
touching distance and at smaller distances. 

There are many nuclei which have positive values of 
the shell-correction energies. Therefore, it is interesting 
to study the effect of the sign of the shell-correction en¬ 
ergy on the nucleus-nucleus potential. The full relaxed- 
density potentials obtained with opposite values =f 4-45 
MeV of the shell-correction energy of 58 Ni and the macro¬ 
scopic part of the potential (the full potential without 
the shell-correction contribution) are presented for the 
system 58 Ni+ 58 Ni in Fig. 4. The shell-correction contri¬ 
bution to the potential is small at large distances and is 
important around the touching point (see Fig. 4). 

According to Eqs. (1) and (15) the depth of the cap¬ 
ture well of the relaxed-density potential with negative 
shell-correction contributions (the case 5E\ + SE 2 < 0) 
becomes shallower than the one of the macroscopic po¬ 
tential. The barrier thickness of the full potential is larger 
than the one for the macroscopic potential. In contrast 
to this, the shell-correction contribution to the potential 
at SEi + 8 E 2 > 0 reduces the barrier thickness and in¬ 
creases the well depth of the full potential in comparison 
to the macroscopic ones. This effect is clearly seen in 
Fig. 4. Because of this effect the deep sub-barrier fusion 
hindrance takes place for a nucleus-nucleus system with 
the strong negative shell-correction contribution to the 
full heavy-ion potential, while the strong positive shell- 
correction contribution to the full potential leads to the 
weak enhancement of the deep sub-barrier fusion cross 
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section; see Ref. [43| for details. 

In conclusion, the phenomenological relaxed-density 
nucleus-nucleus potential with the shell-correction con¬ 
tribution is discussed in detail. The shell-correction con¬ 
tribution to the potential is related to inner structure 
of nuclei, which is disturbed by nucleon-nucleon interac¬ 
tions of colliding nuclei. A simple approach for the evalu¬ 
ation of the shell correction contribution to the full poten¬ 
tial is proposed. The shell-correction contribution shows 
how the full potential for the specific nucleus-nucleus sys¬ 
tem deviates from the global macroscopic potential. The 
global macroscopic potential smoothly depends on A±, 
Z\, A 2 and Zi in contrast to the shell-correction part. 
The shell-correction contribution to the full potential 
is very important at distances smaller than the barrier 


radius. The phenomenological relaxed-density nucleus- 
nucleus potential with the shell-correction contribution 
can reproduce the empirical barrier heights with a value 
of the root mean square error of 0.879 MeV. 
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